COMPUTING OPTIMAL MORSE MATCHINGS 



MICHAEL JOSWIG AND MARC E. PFETSCH 



Abstract. Morse matchings capture the essential structural information of discrete Morse func- 
tions. We show that computing optimal Morse matchings is MV-hard and give an integer pro- 
gramming formulation for the problem. Then we present polyhedral results for the corresponding 
polytope and report on computational results. 



1. Introduction 

Discrete Morse theory was developed by Forman as a combinatorial analog to the classical 

smooth Morse theory. Applications to questions in combinatorial topology and related fields are 
numerous: e.g., Babson et al. 3 , Forman JHI) Batzies and Welker 4 , and Jonsson |2()j . 

It turns out that the topologically relevant information of a discrete Morse function / on a 
simplicial complex can be encoded as a (partial) matching in its Hasse diagram (considered as 
a graph), the Morse matching of /. A matching in the Hasse diagram is Morse if it satisfies 
a certain, entirely combinatorial, acyclicity condition. Unmatched /c-dimensional faces are called 
critical; they correspond to the critical points of rank A; of a smooth Morse function. The total 
number of non-critical faces equals twice the number of edges in the Morse matching. The purpose 
of this paper is to study algorithms which compute maximum Morse matchings of a given finite 
simplicial complex. This is equivalent to finding a Morse matching with as few critical faces as 
possible. 

A Morse matching M can be interpreted as a discrete flow on a simplicial complex A. The flow 
indicates how A can be deformed into a more compact description as a CW complex with one cell 
for each critical face of M. Naturally one is interested in a most compact description, which leads 
to the combinatorial optimization problem described above. This way optimal (or even sufficiently 
good) Morse matchings of A can help to recognize the topological type of a space given as a finite 
simplicial complex. The latter problem is known to be undecidable even for highly structured 
classes of topological spaces, such as smooth 4-manifolds. 

Optimization of discrete Morse matchings has been studied by Lewiner, Lopes, and Tavares |24j . 
Hersh 18 investigated heuristic approaches to the maximum Morse matching problem with ap- 
plications to combinatorics. Morse matchings can also be interpreted as pivoting strategies for 
homology computations; see Furthermore, the set of all Morse matchings of a given simplicial 
complex itself has the structure of a simplicial complex; see • 

Since its beginnings in Lovasz' proof j^S] of the Kneser Conjecture, combinatorial topology 
seeks to solve combinatorial problems with techniques from (primarily algebraic) topology. And, 
conversely, such applications to combinatorics frequently shed light on basic concepts in topology. 
Recent years saw a still growing influx of ideas from differential geometry to the subject. Besides 
discrete Morse theory this concerns, e.g., various notions of discrete curvature which have been 
applied to problems in computational geometry and computer graphics. This is paralleled in the 
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development of a discrete differential geometry. Here the combinatorial cores of known phenomena 
in differential geometry are investigated, while keeping a steady eye on applications in computer 
graphics and mathematical physics; see, e.g., Pinkall and Polthier |2H1, Bobenko and Suris We 
firmly believe that — via discrete Morse theory — techniques from combinatorial optimization are 
highly relevant to these topics. Our contribution is a first step. 

The paper is structured as follows. First we show that computing optimal Morse matchings is 
A/'T^-hard. This issue had been addressed previously by Lewiner, Lopes, and Tavares [^, but we 
believe that their argument has a minor gap. Then we give an integer programming formulation for 
the problem. The formulation consists of two parts: one for the matching conditions and one for 
the acyclicity constraints. This turns out to be related to the acyclic subgraph problem studied by 
Grotschel, Jiinger, and Reinelt jj^l- We derive polyhedral results for the corresponding polytope. 
In particular, we give two different polynomial time algorithms for the separation of the acyclicity 
constraints. The paper closes with computational results. 

Like most of discrete Morse theory, also most of our results extend to arbitrary finite CW- 
complexes. We stick to the simplicial setting, however, to simplify the presentation. 

2. Discrete Morse Functions and Morse Matchings 

We will first introduce discrete Morse functions as developed by Forman Chari ^ showed 
that the essential structure of discrete Morse functions is captured by so-called Morse matchings. 
It turns out that this latter formulation directly leads to a combinatorial optimization problem in 
which one wants to maximize the size of a Morse matching. 

We first need some notation. Let A be a (finite abstract) simplicial complex, i.e., a set of subsets 
of a finite set V with the following property: if F G A and G CI F, then G € A; in other words, A 
is an independence system with ground set V. In the following we will ignore as a member of A. 
The elements in V are called vertices and the sets in A are called faces. The dimension of a face F 
is diuiF := \F\ — 1. Let d = max{dimF : F G 3'} he the dimension of A. We often write z-faces 
for i-dimensional faces. Let 9" be the set of faces of A and let fi = /j(A) be the number of faces 
of dimension i > 0. The maximal faces with respect to inclusion are called facets and 1-faces are 
called edges. The complex A is pure, if all facets have the same dimension. For i^, G € A, we write 
F ~< G F C G and dimF = dimG — 1, i.e., denotes the covering relation in the Boolean 
lattice. The graph of A is the (abstract) graph on V in which two vertices are connected by an 
edge if there exists a 1-face containing both vertices. Throughout this paper we assume that A is 
connected, i.e., its graph is connected. This is no loss of generality since the connected components 
can be treated separately. 

The size of A is defined as the coding length of its face lattice, i.e., if A has n faces, then 
size A = 0{n ■ d ■ logn). Statements about the complexity of algorithms in the subsequent sections 
are always with respect to this notion of size. 

A function / : A — > IR, is a discrete Morse function if for every G E A the sets 

(1) {F:F^G, f{G) < f{F)] and {H:G<H, f{H) < /(G)} 

both have cardinality at most 1. The first set includes the faces covered by face G which are not 
assigned a lower value than G, while the second set includes the faces covering G which are not 
assigned a higher value. The face G is critical if both sets have cardinality 0. A simple example of 
a discrete Morse function can be obtained by setting f{F) = diuiF for every F E A. With respect 
to this function every face is critical. 

Discrete Morse functions are interesting because they can be used to deform a simplicial complex 
into a (smaller) CW-complex that has a cell for each critical face; see Section |31 
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Consider the Hasse diagram H = (3~, A) of A, that is, a directed graph on the faces of A with 
an arc (G, F) E j4 if F -< G; note that the arcs lead from higher to lower dimensional faces. Let 
M C A be a matching in i.e., each face is incident to at most one arc in M. Let H{M) be the 
directed graph obtained from H by reversing the direction of the arcs in M. Then M is a Morse 
matching of A if H{M) does not contain directed cycles, i.e., is acyclic (in the directed sense). 

Chari [Hj observed the following relation between Morse functions and Morse matchings. Let / 
be a discrete Morse function and let M be the set of arcs {G,F) € A such that f{G) < f{F), 
i.e., / is not decreasing on these arcs. A simple proof shows that at most one of the sets in 
can have cardinality one; see Chari [S]. This shows that M is a matching. Since the order given 
by / can be refined to a linear ordering of the faces of A, the directed graph H{M) is in fact 
acyclic and therefore a Morse matching. To construct a discrete Morse function from a Morse 
matching, compute a linear ordering extending H{M) (which is acyclic) and then number the faces 
consecutively in the reverse order. 

Although we loose the concrete numbers attached to the faces when going from a discrete Morse 
function / to the corresponding Morse matching M, we do not loose the information about critical 
faces: Critical faces of / are exactly the unmatched faces of M. Hence, by maximizing \M\ we 
minimize the number of critical faces of /. In fact, the number of critical faces is |3"| — 2 \M\. For 
< j < d, let Cj = Cj{M) be the number of critical faces of dimension j and let c(M) be the total 
number of critical faces. 

It seems helpful to briefly describe the case of Morse matchings for a one-dimensional simplicial 
complex A. Then A represents the incidences of a graph G. A Morse matching M of A matches 
edges with nodes of G. Let G be the following oriented subgraph of G: take all edges which are 
matched in M and orient them towards its matched node. Since M is a matching, this construction 
is well defined and the in-degree of each node is one. The acyclicity property shows that G contains 
no directed cycles and hence is a branching, i.e., the underlying graph is a forest and each (weakly) 
connected component has a unique root. Therefore, the Morse matchings on a graph G are in 
one-to-one correspondence with orientations of subgraphs of G which are branchings. 

Generalizing this idea, Lewiner, Lopes, and Tavares developed a heuristic for computing 
optimal Morse matchings. This heuristic computes maximum Morse matchings for combinato- 
rial 2-manifolds, i.e., Morse matchings with maximal cardinality. However, for general simplicial 
complexes this problem is A/'P-hard, see Section |^ 

3. Properties of Morse Matchings 

In this section we briefly review some important properties of Morse matching which we need in 
the sequel. 

Let F be a facet of A and let G be a facet of F, which is not contained in any other facet of A. 
The operation of transforming A to A \ {F, G} is called a simplicial or elementary collapse. We 
will simply use collapse in the following. 

Proposition 3.1 (Forman 0). Let A be a simplicial complex and S a subcomplex of A. Then 
there exists a sequence of collapses from A to T, if and only if there exists a discrete Morse function 
such that A \ S contains no critical face. 

Forman (Qj also proved the following result, which describes one of the most interesting features 
of Morse matchings: 

Theorem 3.2. Let A be a simplicial complex and M be a Morse matching on A. Then A is 
homotopy equivalent to a CW- complex containing a cell of dimension i for each critical face of 
dimension i. 
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We refer to Munkres [22| for more information on CW-complexes. By this result we can hope for 
a compact representation of the topofogy of A (up to homotopy) by computing a Morse matching 
with few critical faces. This is the main motivation for the combinatorial optimization problem 
studied in this paper. 

Let K he a field and let f3j = f3j{K) to be the Betti number for dimension j over K for A; see 
again Munkres [23 for details. Forman 9 proved the following bounds on the number of critical 
faces Cj of a Morse matching M: 

Theorem 3.3 {Weak Morse inequalities). Let K he a field, A be a simplicial complex, and M a 
Morse matching for A. We have 

(2) Cj > Pj for allj = 0,...,d 
and 

(3) CO - Ci + C2 - • • • + i-lfcd = /?o - A + /?2 - • • • + i-iyPd- 

The Betti numbers over Q and finite fields can easily be obtained in polynomial time (in the size 
of A), by computing the ranks of the boundary matrices for each dimension. Although harder to 
compute (see Iliopoulos |19j). the homology over Z can be used to choose among the finite fields 
or Q, in order to obtain the strongest form of the Morse inequalities 

4. Hardness of Optimal Morse Matchings 

In this section we prove TVP-hardness of the problem to compute a maximum Morse matching, 
i.e., to find a Morse matching M with maximal cardinality. As we saw previously, this is equivalent 
to minimize the number of critical faces. 

We want to reduce the following collapsihility problem, introduced by Egecioglu and Gonzalez |Sj, 
to the problem of finding an optimal Morse matching: Given a connected pure 2-dimensional 
simplicial complex A, which is embeddable in I^? and an integer k, decide whether there exists 
a subset K, of the facets of A with |/C| < k such that there exists a sequence of collapses which 
transforms A\/C to a 1-dimensional complex. Egecioglu and Gonzalez proved that this collapsibility 
problem is strongly AAP-complete. Using Proposition 13.11 this results reads in terms of discrete 
Morse theory: 

Theorem 4.1. Given a connected pure 2-dimensional simplicial complex A, which is embeddable 
in I^? , and a nonnegative integer k, it is MV- complete in the strong sense to decide whether there 
exists a Morse matching with at most k critical 2-faces. 

When k is fixed, we can try all possible sets /C of size at most k and then decide whether the 
resulting complex is collapsible to a 1-dimensional complex in polynomial time. Therefore we let k 
be part of the input. 

We need the following construction. Consider a Morse matching M for a simplicial complex A, 
with dim A > 1. Let T{M) be the graph obtained from the graph of A by removing all edges 
(1-faces) matched with 2-faces. Note that T{M) contains all vertices of A. 

Lemma 4.2. The graph T(AI) is connected. 

Proof. Without loss of generality we assume that dim A > 2. Otherwise T{M) coincides with the 
graph of A, which is connected (recall that A is connected). 

Suppose that T(M) is disconnected. Let be the set of nodes in a connected component of 
r(M), and let C be the set of cut edges, that is, edges of A with one vertex in A^ and one vertex 
in its complement. Since A is connected, C is not empty. By definition of T(M), each edge in C is 
matched to a unique 2-face. 
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Consider the directed subgraph D of the Hasse diagram consisting of the edges in C and their 
matching 2-faces. The standard direction of arcs in the Hasse diagram (from the higher to the 
lower dimensional faces) is reversed for each matching pair of M, i.e., D is a subgraph of H{M). 

We construct a directed path in D as follows. Start with any node of D corresponding to a cut 
edge ei. Go to the node of D determined by the unique 2- face ri to which ei is matched to. Then ri 
contains at least one other cut edge 62, otherwise ei cannot be a cut edge. Now iteratively go to 62, 
then to its unique matching 2-face r2, choose another cut edge 63, and so on. We observe that we 
obtain a directed path ei, ri, 62, T2, . . . in L), i.e., the arcs are directed in the correct direction. 

Since we have a finite graph at some point the path must arrive at a node of D which we 
have visited already. Hence, D (and therefore also H{M)) contains a directed cycle, which is a 
contradiction since M is a Morse matching. □ 

Now pick an arbitrary node r and any spanning tree of T{M) and direct all edges away from r. 
This yields a maximum Morse matching on r(M); see end of Section |2J It is easy to see that 
replacing the part of M on T{M) with this matching yields a Morse matching. This Morse matching 
has only one critical vertex (the root r). Note that every Morse matching contains at least one 
critical vertex; this can be seen from the Morse inequalities Theorem 13.31 or directly by observing 
that we get a directed cycle in the Hasse diagram if every vertex is matched to an edge. Furthermore, 
the total number of critical faces can only decrease, since we computed an optimal Morse matching 
on r(M). The number of critical i-faces for i > 2 stays the same. We have thus proved the 
following, which is also implicit in For man P|. 

Corollary 4.3. Let M he a Morse matching on A. Then we can compute a Morse matching M' 
in polynomial time which has exactly one critical vertex and the same number of critical faces of 
dimension 2 or higher as M , such that c{M') < c(M). 

We can now prove the hardness result. 

Theorem 4.4. Given a simplicial complex A and a nonnegative integer c, it is strongly J\fV- 
complete to decide whether there exists a Morse matching with at most c critical faces, even if A 
is connected, pure, 2 -dimensional, and can be embedded in R^. 

Proof. Clearly this problem is in AfV. So let (A, k) be an input for the collapsibility problem. 
We claim that there exists a Morse matching with at most k critical 2-faces if and only if there 
exists a Morse matching with at most g{k) := 2{k + 1) — x(^) critical faces altogether. Here, 
x(A) = /?o — /?! + • • • + {—l)'^Pd is the Euler characteristic, which can be computed in polynomial 
time; see Section |31 Hence g is a polynomial-time computable function. Using Theorem 14. II then 
finishes the proof. 



So assume that M is a Morse matching on A with at most k critical 2-faces. We use Corollarv 14.31 
to compute a Morse matching M', in polynomial time, such that co(M') = 1, C2(M') = C2{M), 
and c(M') < c(M). By the Morse equation of Theorem lOl we have ci(M') = C2(M') + 1- x(A). 
Since c(M') = co(M') + ci(M') + C2(M') it follows that 



Solving for c(M'), it follows that M' has at most 2{k + 1) — x(^) critical faces altogether. 

Conversely, assume that there exists a Morse matching M with at most g{k) critical faces. 
Computing M' as above, we obtain by (HJ, that 



(4) 



C2(M) = C2(M') = i(c(M') + x(A))) - 1. 



C2(M) =C2(M') < i 



l{g{k)+x{A)))-l = k, 



which completes the proof. 



□ 
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Lewiner, Lopes, and Tavares [2^ showed that it is A/''P-hard to compute an optimal Morse 
matching with exactly one critical vertex. They claimed it for the general case, but we do not see 
an argument similar to Lemma l4.2l above. 

Since there exists a Morse matching with at most c critical faces if and only if there exists a 
Morse matching of size at least — c), we proved: 

Corollary 4.5. Let A be as in Theorem \4-4\ ^''T-d m be a nonnegative integer. Then it is AfV- 
complete in the strong sense to decide whether there exists a Morse matching of size at least m. 

We do not know about the complexity status for this problem with m fixed. 

Egecioglu and Gonzalez even proved that the collapsibility problem is as hard to approximate 
as the set covering problem. In particular, the collapsibility problem cannot be approximated better 
than within a logarithmic factor in polynomial time, unless V = MV. Using this, Lewiner, Lopes, 
and Tavares |24j claimed that the problem to compute a Morse matching minimizing the number 
of critical faces is hard to approximate. However, the function g used in the proof above is not 
"approximation preserving" and we do not see how the non-approximability result carries over. 

Similarly, the problem to approximate the size of a Morse matching seems to be open. 

5. An IP-Formulation 

In this section we introduce an integer programming formulation for the problem to compute a 
Morse matching of maximal size. It can easily be extended to arbitrary weights. 

We use the following notation. We depict vectors in bold font. Let Cj be the zth unit vector and 
let 1 be the vector of all ones. For any vector x S IR" and / C {1, . . . , n} we define 

Furthermore, for S C {1, . . . n}, I{S) E K," denotes the incidence vector of S. 

For a node f in a directed graph, let 6{v) be the arcs incident to v, i.e., the arcs having v as one 
of their endnodes. For a subset A' C A, we denote by A^(^') the nodes incident to at least one 
arc in A' . Throughout this article, all directed or undirected cycles are assumed to be simple., i.e., 
without node repetitions. 

For ease of notation, we consider the Hasse diagram H as directed or undirected depending on 
the context; we will explicitly use directed when we refer to the directed version. 

We split H into d levels Hq = {3^,Aq), . . . , Hd-i = {S^'^'^ , Ad-i) , where Hi denotes the level 
of the Hasse diagram between faces of dimension i and i + 1. Then A is the disjoint union of 
• • • ) Ad-i and 3^~^ n 3^ consists of the faces of dimension i. Recall that the arcs in the Hasse 
diagram are directed from the higher to the lower dimensional faces. 

Let M C j4 be a Morse matching of A. By definition, its incidence vector x = I{M) G {0, 1}"^ 
satisfies the matching inequalities: 

(5) x{5{F))<l VFgJ. 

Now assume that for some M ^ A there exists a directed cycle D in H{M). Then in D "up" 
and "down" arcs alternate; see for example Figure E In particular, the size of D is always even. 
Hence, \\D\ arcs are contained in M, i.e., are reversed in H{M). Furthermore we note: 

Observation. Let M C ^ be a matching. If Z) is a directed cycle in H{M), the edges in D can 
only belong to one level Hi (i G {0, . . . ,d - 1}), i.e., {dimF : F G N{D)} = + 1}. 

Putting these arguments together we obtain: If M is acyclic, x = I{M) satisfies the following cycle 
inequalities: 

(6) x{C)<\\C\-l VCGCi, i = l,...,d-l. 
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Figure 1: Example for a directed cycle of size 6; at least three arcs with reversed orientation (pointing "up") are 
necessary to close a 6-cycle in the Hasse diagram of a simplicial complex. 



where Ci are the cycles in Hi. 

Conversely, it is easy to see that every x G {0, 1}"^ which fulfills inequalities © and © is the 
incidence vector of a Morse matching. Hence, we arrive at the following IP formulation for the 
problem to find a maximum Morse matching: 

(MaxMM) max l^x 

s.t. x{6{F)) < 1 VF G J 

x{C)<^\C\-l \/ C £Ci, i = l,...,d-l 

X G {0,1}^. 

We define the corresponding polytope: 

Pm = conv {a; G {0, 1}^ : x satisfies © and (HI)}. 

Let M be a Morse matching and x = I{M) be its incidence vector. Then F G 9" is a critical face 
with respect to M if and only if it is unmatched by M, i.e., x(6{F)) = 0. Hence, the total number 
of critical faces is: 

(7) c(M) = ^(i- xa) = m- 2Y,xa = m-2i^x, 

since every arc is incident to exactly two nodes. Using this formula one can easily switch between 
the number of critical faces and the number of arcs in a Morse matching. 

The LP relaxation of MaxMM can be strengthened by using the weak Morse inequalities ((2) of 
Theorem 13.31 Applying ((T)), this yields the following Betti inequality for dimension i: 

(8) Yl E ^-)^f^^ ^ EE ^a<f^-P^. 
F:dimF=i ae5(F) F:dimF=i a&S{F) 

Observe that we can choose the field in Theorem 13.31 to employ the Morse inequalities in their 
strongest form. 

Note. The cycle inequalities ^ are similar to the cycle inequalities for the acyclic subgraph prob- 
lem (ASP); see Jiinger [22j, and Grotschel, Jiinger, and Reinelt JHI- The separation problem 
for however, is more complicated than the corresponding problem for ASP; see Section [5.21 

Furthermore, there is a similarity to the relation between the ASP and the linear ordering 
problem (see Reinelt |29j . and Grotschel, Jiinger, and Reinelt jl4j): an alternative formulation for 
our problem can be obtained by adding matching inequalities to the linear ordering formulation; 
this directly models discrete Morse functions as linear orderings of the faces. Since this formulation 
is based on the relation between faces, it leads to quadratically many variables in the number 
of faces; therefore we have opted for the above formulation, at the cost of having to solve the 
separation problem for the cycle inequalities; see Section EH 
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Figure 2: Example for a non-monotone behavior of acyclic matchings. The directed graph on the right, obtained 
from the left graph by reversing the dashed arcs, is acyclic. However, if the top arc is set to its original orientation, 
the graph is not acyclic anymore. This shows that subsets of acyclic matchings are not necessarily acyclic. 



5.1. Facial Structure of Pm- It is easy to see that Pm is a full dimensional polytope and 
Xa ^ defines a facet for every a £ A. Furthermore, Pm is monotone, since every subset of a 
Morse matching is a Morse matching. It is well known that this implies that every facet defining 
inequality ot^x < (i not equivalent to the non-negativity inequalities fulfills: a > 0, /? > 0; see 
Hammer, Johnson, and Peled [12|. 

Interestingly, if we generalize Morse matchings to acyclic matchings for arbitrary graphs, the 
collection of such acyclic matchings is not necessarily monotone anymore; see the example in 
Figure 121 Therefore, the structure of the generalized problem is likely to be more complicated. 

We have the following two results: 

Proposition 5.1. The matching inequalities x{6{F)) < 1 define facets of Pm for F £ T, except if 
\5{F) \ = 1, i.e., F is a vertex. 

Proof. Let F be a face with > 1 (note that = does not occur). We can assume 

that A = {ai, . . . , a^, flfc+i, • • • , am}, where 5{F) = {oi, . . . , a^}. For i = A; + 1, . . . , m, observe 
that ai cannot be adjacent to every arc in 5{F): since > 1, would either be incident to at 

least two nodes of the same dimension or to two nodes whose dimensions are two apart, which is 
impossible. Therefore, choose p{i) G {1, . . . , fc} such that and ap(j) are not adjacent. It follows 
that ej + ep(j) G Pm ■ Then 

) • • • 1 ^fc ) ^fc + l ~l" ^p(fc-l-l) ) • • • ; ~l" ^p(m) 

are affinely independent and fulfill x{5{F)) = 1. □ 

It follows that the inequalities < 1, a € A, never define facets. 
Theorem 5.2. The cycle inequalities © define facets of Pm- 

Proof. We extend the corresponding proof by Jiinger for the ASP. 

Let C be a cycle in H. Without loss of generality assume that A = {ai, . . . , a^, a/j+i, . . . , am}, 
where C = {ai, . . . ,ak) and k is even. We will construct affinely independent feasible vectors 
1?!, . . . , Ufc+i, . . . , Vm satisfying the cycle inequality corresponding to C with equality. 

Let Ci = {ai, 03, . . . , Ok-i} and C2 = {02, 04, . . . , a^}. Hence Ci and C2 are the "up" and 
"down" arcs in C. 

Define 

\HC2\la,}) itajeCz 

Hence, for i = 1, . . . , /c we have Vi{C) = | — 1. 

For i = k + 1, . . . ,m, consider Oj = {u, v} ^ C. We have four cases. 
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Figure 3: Illustration of the first case in the proof of Theorem 15.21 The sets Pi and P2 are shown by continuous 
lines. The edges in Ci are drawn gray and hence Pi C Ci; edges in C2 are drawn black. The dashed edges incident 
to u and v are not considered. The right hand side shows the graph embedded in the Hasse diagram. 



> u,v £ N{C): Let C := C\ [5{u)u5{v)) . We have that \C\ = k — A (since there exist no odd cycles) 

and C spUts into two odd nonempty parts C\ and C2, which are both paths. Let ki := \Ci \ and 
^2 := \C2\', ki and k2 are odd, since u and v are on opposite sides of the bipartition. We choose 
a subset Pi C Ci by taking every second arc in order to get |-Pi| = similarly we choose 

P2 C C2 with IP2I = -^^^ By construction either Pi C Ci or Pi C C2 and either Pj n C2 = 
or Pj n Ci = for z = 1, 2. An easy calculation shows that |Pi U P2I = I — 1. See Figure El for 
an illustration of this case. Then define Vi := /(Pi U P2 U {oj}). 

> u ^ C, V £ C: Here we define Vi := I{Ci \ 6{v) U {oj}). 
n e C, V ^C: Define Vi := I{Ci \ 5{u) U {aj). 

> u,v ^ C: Choose any a € Ci and define Vi := I{Ci \ {a\ U {fli}). 

It is easy to check in each case that Vi G Pm and that Vi{C) = -I — 1. 

It can be shown that the m vectors vi,. . . ,Vm are affinely independent, which concludes the 
proof. □ 

The separation problem for the cycle inequalities is discussed in the next section. 

5.2. Separating the Cycle Inequalities. Of course, there are exponentially many cycle inequal- 
ities (jnj. Hence we have to deal with the separation problem for these inequalities. 

For the separation problem, we can assume that we are given x* € [0,1]"^, which satisfies all 
matching inequalities ©• We consider the separation for each graph Hi in turn, i = 0, . . . , d — 1. 
The problem is to find an undirected cycle C in Hi such that 

x*{C) > i|Cl -1 

or conclude that no such cycle exists. In the next sections we describe two methods to solve this 
separation problem in polynomial time. 

5.2.1. Undirected Shortest Path with Conservative Weights. A usual trick to solve the above sep- 
aration problem is to apply an affine transformation and obtain a shortest cycle problem. The 
transformation suitable for our needs is a;' = ^1 — a;, which yields: 

x{C)<^\C\-l ^ x'{C)>l. 

The separation problem can now be solved as follows: compute a shortest cycle in Hi with respect 
to the weights ^1 — x*. If its weight is at most 1, this cycle yields a violated cycle inequality, 
otherwise no such cycle exists. 

However, the weights can be negative and we have to rule out negative cycles in order to apply 
polynomial time methods from the literature; that is, we want the weights to be conservative. 
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W2 

O 



us 



W3 

o 



U4 




({U2,U4}, W2) 

{{u-i,U4,},Wz) 



Figure 4: Example of the construction in Section fS. 2. 21 Left: original graph G. Right: constructed graph G' . The 
6-cycle on the left corresponds to the 3-cycle on the right (both shown with dashed lines). 



Lemma 5.3. There exists no cycle of negative weight in Hi with respect to t^^—x*, forO <i< d—1. 

Proof. Let C = (ai, . . . , a^) be a cycle in Hi and let Fi, . . . ,Fk be the faces that are visited by C. 
Recall the x* satisfies the matching inequalities. We obtain 

(9) E E < = = 2x*(C), 

j=la€5{F,)nC a€C 

since each edge weight is counted twice in the first term. Applying the Matching inequalities © on 
the left hand side yields that x*{C) < = ^\C\. Hence, the weight of C with respect to ^1 — x* 
can be bounded as follows 

Y^{i-x:) = i\c\-x*{c)>o, 

which proves the lemma. □ 

We have now reduced the separation problem to finding a shortest cycle in a weighted undirected 
graph G = {V, E) without negative cycles. 

By using T-join techniques, one can compute a shortest path in an undirected graph without neg- 
ative cycles in 0{ni{mi + Ui logrij)) time, where in this formula Ui = and rrii = \ Ai\] see Schri- 
jver Chapter 29]. It follows that a shortest cycle can be computed in 0{mini{mi + njlognj)) 
time. Since \Ai\ < (i + 2)nj, this leads to an 0(^{d + l)'^n^ + (d + l)n'^logn) overall algorithm, 
where n := I?'! is the number of faces and d is the dimension of the complex. 

5.2.2. Transforming the Graph. Another method for the separation problem of cycle inequalities, 
which is easier to implement, works as follows. 

Let G = {U (JW, E) he a. bipartite graph, e.g., G = Hi (z € {0, . . . , d — 1}), the z-th level of the 
Hasse diagram. Let i : E ^ lR.>o be a length function for the edges of G. In the following we write 
i{u, v) = £{v, u) for the length £{{u, v}). 

We construct a graph G' = (y',E') and lengths i' : E' ^ lR.>o ^ follows; see Figure |1] for an 
example. The set of nodes of G' is 

{{{u,u'},w) : {u,u'} CU, weW, {u,w} G E, {u',w} G E}. 

We have an edge between two nodes {{ui,u'i},wi) and ({n2, 1*2}, 11)2) if 

{ui, u'l} n {u2, U2} 7^ and wi 7^ ^2. 

The length of such an edge e' is defined by 

e'{e') = \ {i{ux,w{) + ^K, + £{U2,W2) + ^K, «;2)) . 

We now consider the relation of cycles in G and G' . 
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Lemma 5.4. C = {uq^wojUijIVi, . . . ,Wk-i,ui) is a cycle in G with k > 1 of length i{C) if and 
only if 

C = {{{uo,Ui},Wo), {{ui,U2},Wi), {{Uk-l,Ui},Wk-l), {{uo,ui},wo)) 

is a cycle in G' with f (C) = i[C). 

Proof. In the following we compute with indices modulo k. 

First observe that C' is well defined: Each ({wj, ttj+i}, Wj) is a node of G' , since {ui,Wi}, 
{uiJ^i^Wi} E E. Furthermore, 

{{{ui,Ui+i},Wi), {{ui+i,Ui+2},Wi+i)] € E' 

since Wi ^ Wi+i and {nj,Mj+i} n {uj+i, 114+2} 7^ (because k > 1). It is a cycle since k>2. 
The weight of C' can be calculated as follows: 



k 
k 

= {i{ui,Wi)+£{ui+i,Wi)) =£{C), 



i=0 

where the second to last equation follows since every edge in C occurs twice in the summation. □ 

The previous lemma does not cover cycles in G of length four. These do not occur for the case 
of G = Hi, since Hi is a level in the Hasse diagram of a simplicial complex. Moreover, cycles of 
length four can readily be detected in the construction of G' and handled accordingly (there is only 
a polynomial number of them). 

To solve our separation problem, let G = Hi, i € {0, . . . ,d — 1}, and £{e) = x* for e G G. 
Then we have i'{e') € [0, 1] for each e' G E' , because of the matching inequalities. We now set 
i{e') = 1 — i'{e') for e' G G' and hence i{e') G [0, 1]. Let C be a cycle in G with at least six edges 
and C' be the corresponding cycle in G'. Note that \G'\ = ^\G\. We have: 



iiC')= i{e')= J](l-/(e'))<l 

e'GC e'GC 

^ Y ^'(e') > \C'\ - 1 



e'eC 

^ e'{c') > \c'\ - 1 

^ i{C)>l\C\-l (by Lemma E31). 

Hence, C violates the cycle inequality © if and only if £{C') < 1. Since £{e') > 0, we can use 
the Floyd- Warshall algorithm to solve the separation problem in time 0(|yp); see Korte and 
Vygen |^. 

liG = Hi and W is the part arising from the higher dimensional faces, we have \V'\ = {'\^)\W\ = 
This leads to an ©((d -|- l)^n^) algorithm for separating cycle inequalities, which is 
roughly as fast as the method discussed in Section 15.2.11 but much easier to implement. 
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Table 1: Computational results of the branch-and-cut algorithm with separating cycle inequalities and Gomory cuts. 



name 


n 


ni 


d 


nodes 


depth 


time 





c 


solid_2_torus 


24 


42 


2 


1 





0.00 


2 


2 


simon2 


31 


60 


2 


1 





0.00 


1 


1 


projective 


31 


60 


2 


1 





0.01 


3 


3 


bjorner 


32 


63 


2 


1 





0.05 


2 


2 


nonextend 


39 


77 


2 


6 


5 


0.16 


1 


1 


simon 


41 


82 


2 


1 





0.18 


1 


1 


dunce 


49 


99 


2 


385 


10 


2.62 


1 


3 


c-ns3 


63 


128 


2 


349 


10 


3.47 


1 


3 


c-ns 


75 


152 


2 


28 


10 


1.95 


1 


3 


c-ns2 






z 


1 A 
14 


7 


111 
1.11 


1 


1 


zieffler 


119 


310 


3 


1 





0.01 


1 


1 


QTi 1 pnhan in 


167 


434 


3 


1 





25.24 


1 


1 


lockeberg 


216 


600 


3 


1 





36.25 


2 


2 


rudin 


215 


578 


3 


77 


30 


103.78 


1 


1 


mani-walkup-D 


392 


1112 


3 


111 


23 


512.81 


2 


2 


mani-walkup-C 


464 


1312 


3 


135 


83 


1658.02 


2 


2 


MNSB 


103 


267 


3 


12 


10 


73.39 


1 


1 


MNSS 


250 


698 


3 


292 


110 


750.36 


2 


2 


CP2 


255 


864 


4 


230 


80 


558.14 


3 


3 



6. Computational Results 

In this section we report on computational experience with a branch-and-cut algorithm along 
the lines of Section |SI The C-I-+ implementation uses the framework SCIP (Solving Constraint 
Integer Programs) by Achterberg, see |lj. It furthermore builds on po lymake; see HSHlH]. As an 
LP solver we used CPLEX 9.0. 

As the basis of our implementation we take the formulation of MaxMM in Sectional Matching 
inequalities © and Betti inequalities (jEJ (together with variable bounds) form the initial LP. Cycle 
inequalities © are separated as described in Section 15.2.21 Additionally, Gomory cuts are added. 
As a branching rule we use reliability branching implemented in SCIP, a variable branching rule 
introduced by Achterberg, Koch, and Martin [2]. 

We implemented the following primal heuristic. First a simple greedy algorithm is run: We 
start with the empty matching M = 0. We add arcs of the Hasse diagram to M in the order of 
decreasing value of the current LP solution as long as M stays an acyclic matching (which can 
easily be tested). Then the outcome is iteratively improved by a method described in Forman j^]: 
One searches for a unique path between two critical faces in H{M). Such a path is alternating 
with respect to M. Then M can be augmented along the path (the new matching is the symmetric 
difference of M and the path). As is easily seen, this generates an acyclic matching, because the 
path is unique. This heuristic turns out to be extremely successful; see below. 

We tested the implementation on a set of simplicial complexes collected by Hachimori; see for 
more details. Additionally, we considered the following complexes: CP2 (complex projective plane), 
CP2+CP2 (connected sum of CP2 with itself), MNSB and MNSS ((vertex) minimal non-shellable bah 
and sphere, respectively; see Lutz i26i ). 

All computational experiments were run on a 3 GHz Pentium machine running Linux. In the 
tables of computational results, n denotes the number of faces, m the number of arcs in the Hasse 
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Table 2: Computational results of the branch-and-cut algorithm without separation. 



name 


n 


ni 


d 


nodes 


depth 


time 




H 


c 


solid_2_torus 


24 


42 


2 


1 





0.00 


2 


2 


simon2 


31 


60 


2 


1 





0.01 


1 


1 


projective 


31 


60 


2 


1 





0.00 


3 


3 


bjorner 


32 


63 


2 


1 





0.01 


2 


2 


nonextend 


39 


77 


2 


3 


2 


0.02 


1 


1 


simon 


41 


82 


2 


4 


3 


0.02 


1 


1 


dunce 


49 


99 


2 


168367 


42 


145.60 


1 


3 


c-ns3 


63 


128 


2 


3665581 


53 


3940.40 


1 


3 


c-ns 


75 


152 


2 


16625713 


58 


19359.69 


1 


3 


c-ns2 


It) 




o 
z 


A 

4 


Q 

o 


n HQ 
U.Uo 


i 


1 


zieffler 


119 


310 


3 


1 





0.01 


1 


1 


QTi 1 pnhan in 


167 


434 


3 


21 


20 


0.68 


1 


1 




216 


600 


3 


1 





0.05 


2 


2 


rudin 


215 


578 


3 


81 


80 


3.18 


1 


1 


mani-walkup-D 


392 


1112 


3 


107 


100 


2.00 


2 


2 


mani-walkup-C 


464 


1312 


3 


1498 


456 


30.54 


2 


2 


MNSB 


103 


267 


3 


1 





0.01 


1 


1 


MNSS 


250 


698 


3 


163 


126 


4.63 


2 


2 


CP2 


255 


864 


4 


198 


190 


4.77 


3 


3 


CP2+CP2 


460 


1592 


4 


5178 


534 


110.21 


4 


4 



diagram (= number of variables), d the dimension, nodes the number of nodes in the branch-and- 
bound tree, depth the maximal depth in the tree, time the computation time in seconds, (3 is the 
lower bound obtained by adding all Betti inequalities 0, and c the number of critical faces in the 
optimal solution. 

Our implementation could not solve the larger problems of Hachimori's collection in reasonable 
time: bing, knot, poincare, nonpl_sphere, and nc_sphere. In fact, for poincare we ran our code 
in different settings, each for about a week - without success. 

Table ^ shows the results of a computation where we separate cycle inequalities and Gomory 
cuts and run the heuristic every 10th level. At most seven separation rounds of cycle inequalities 
were performed at a node. We do not report results on the problems by Moriyama and Takeuchi in 
Hachimori's collection - they all could be solved within a second. The version with cut separation 
could not solve CP2+CP2 within 90 minutes. 

For most problems the bound obtained by adding Betti inequalities (|H1), as indicated in col- 
umn "/?", is tight. This means that the algorithm is done once an optimal solution is found. This 
usually happens very fast and shows that the heuristic is efficient. In fact, there are only three 
problems for which the bound is not tight and could be solved by our algorithm (dunce, c-ns, and 
c-ns3). These three problems are solved easily by the version with cut separation. In our problem 
set there exists no hard but still solvable problem with a "Betti bound" which is not sharp. We 
can therefore not estimate the limits of our implementation for these cases (poincare is the next 
larger problem of this kind with 1112 variables, but we could not solve it). 

The tractability of problems with a tight "Betti bound" is supported by the results obtained 
by running the implementation without any separation; see Table |5J Only integer solutions are 
checked whether they are acyclic and the heuristic is run every 10th level. This essentially is a test 
of the performance of the primal heuristic. Indeed, all problems with tight "Betti bound" were 
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solved within a few seconds (CP2+CP2 and mani-walkup-C being the exception, but could be solved 
within two minutes). The results for the problems c-ns, c-ns3, and dunce show that the cycle 
inequalities and Gomory cuts are very effective in reducing the number of nodes in the tree and 
the computing time for problems where the "Betti bound" is not sharp. 

Summarizing, we can say that our implementation can solve large instances with up to about 
1500 variables if the bounds from the Betti numbers are tight and small instances with up to about 
150 variables if the bounds are not tight. In all the instances computed so far, the topology of the 
spaces involved was known. In the future, we plan to apply our techniques to other cases. 

Acknowledgment. We are indebted to Tobias Achterberg for his support of the implementation. 
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